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1 Introduction 



The last decade of quark flavour experiments has shown an impressive success of the simple CKM 
mechanism for flavour mixing and CP violation: All measurements of rare decays (AF = 1), of mixing 
phenomena (AF = 2), and of all CP violating observables at tree and loop level have been consistent 
with the Cabibbo-Kobayashi-Maskawa (CKM) theory of the Standard Model (SM); in other words 
none of the former and present flavour experiments including the first generation of the B factories 
at KEK (Belle experiment at the KEKB e + e~ collider) fl] and at SLAC (BaBar experiment at the 
PEP-II e + e~ collider) [2j, and the Tevatron B physics programs (CDF [3] and DO [4| experiments) 
has found an unambiguous sign of New Physics (NP). Moreover, all first results based on the high 
statistics of the LHCb experiment |5] are again very well in agreement with the CKM theory of the 
SM. Of course there have been and there are still so-called tensions, anomalies, or puzzles in the quark 
flavour data at l-,2-, or 3-<r level, however, until now they all have disappeared after some time when 
more statistics had been collected. 

This means that flavour- violating processes between quarks are governed by a 3 x 3 unitarity matrix 
referred to as the CKM matrix [6j[7]. In particular, the one phase among the four real independent 
parameters of the CKM matrix represents the dominating source of CP violation and it allows for a 
unified description of all the CP violating phenomena in the SM. This success of the CKM theory of 
CP violation was honored with the Nobel Prize in Physics in 2008. 

This feature is somehow unexpected because in principle (loop-induced) flavour changing neutral 
current (FCNC) processes like B — > X s j offer high sensitivity to NP; additional contributions to the 
decay rate, in which SM particles are replaced by new particles such as the supersymmetric charginos 
or gluinos, are not suppressed by the factor a/47r relative to the SM contribution. Thus, FCNC 
decays provide information about the SM and its extensions via virtual effects to scales presently 
not accessible (for reviews see Refs. [8j[9] ) . This is complementary to the direct production of new 
particles at collider experiments 



Within this indirect search for NP there is a an ambiguity of the NP scale. In the model- 
independent approach using the effective electroweak hamiltonian, the contribution to one specific 
operator Oi can be parametrized via (Cg M / Mw + C* NP / Anp) x Oi where the first term represents 
the SM contribution at the electroweak scale M\y and the second one the NP contribution with 
an unknown coupling Cg M and an unknown NP scale Anp- The non-existence of large NP effects in 
FCNC observables in general implies the infamous flavour problem, namely why FCNC are suppressed: 
Either the mass scale of the new degrees of freedom Anp is very high or the new flavour-violating 
couplings C NP are small for (symmetry?) reasons that remain to be found. For example, assuming 
generic new flavour-violating couplings of 0(1), the present data on K-K mixing implies a very high 
NP scale of order 10 3 -10 4 TeV depending on whether the new contributions enter at loop- or at tree- 
level. In contrast, theoretical considerations on the Higgs sector, which is responsible for the mass 
generation of the fundamental particles in the SM, call for NP at order 1 TeV. As a consequence, any 
NP below the 1 TeV scale must have a non-generic flavour structure. 

The hypothesis of minimal flavour violation (MFV) [12 - 14 , is a formal model-independent solution 
to the NP flavour problem. It assumes that the flavour and the CP symmetry are broken as in the 
SM. Thus, it requires that all flavour- and CP-violating interactions be linked to the known structure 
of Yukawa couplings. A renormalization-group invariant definition of MFV based on a symmetry 
principle is given in Ref. [14] ; this is mandatory for a consistent effective field theoretical analysis of 
NP effects. 

The MFV hypothesis is far from being verified. There is still room for sizable new effects, and new 
flavour structures beyond the Yukawa couplings are still compatible with the present data because 
the flavour sector has been tested only at the 10% level especially in the b — > s transitions. However, 
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the MFV hypothesis represents an important benchmark in the sense that any measurement which is 
inconsistent with the general constraints and relations induced by the MFV hypothesis unambiguously 
indicates the existence of new flavour structures. 

This implies the main purpose of the present paper, namely to derive the consequences of the 
MFV hypothesis based on the latest LHCb data. This was done some time ago in Ref. |15| for 
AF = 1 observables. In particular the impressive start of the LHCb experiment suggests to update 
this analysis. 

Besides the new data from the B factories and from the LHCb experiment we implement some 
additional theoretical improvements compared to the previous analysis of AF = 1 processes: The 
exclusive B — > K*£ + £~ decay was analyzed within a simple form factor analysis. Here we use the up- 



to-date theoretical tools of QCD-improved factorization for the low-t/ 2 region 16, 17 and the recently 
proposed OPE methods for the high-q 2 region 18,19. Secondly, we skip the approximation that 



the NP contributions to the electromagnetic and chromomagnetic operators appear in a fixed linear 
combination which was necessary in the previous analysis due to the limited number of independent 
experimental measurements. 

The MFV analysis of AF = 2 mixing phenomena decouples from the AF = 1 analysis within the 
standard MFV framework. The AF = 2 analysis in Ref. [20] was recently updated in Refs. |21[|22| . 

The paper is organized as follows: In Section 2 we recall the definition of the RG invariant 
definition of the MFV hypothesis and in particular the effective hamiltonian within this framework. 
In Section 3 we work out the dependence on the non-standard Wilson coefficients of the MFV effective 
theory for all the flavour observables used in our analysis. We also discuss the various sources of the 
uncertainties in the theoretical predictions. In Section 4 we give some numerical details and in Section 
5 we discuss our results. 



2 Effective hamiltonian with MFV 

2.1 MFV hypothesis 

The SM gauge interactions are universal in quark flavour space, this means the gauge sector of the 
SM is invariant under the flavour group Gfl avouI . which can be decomposed as 

G flavour = U(3) Ql x U{3) Ur x U(3) Dr . (2.1) 

In the SM this symmetry is only broken by the Yukawa couplings. Any new physics model in 
which all flavour- and CP-violating interactions can be linked to the known Yukawa couplings is 
minimal flavour violating. In order to implement this principle in a renormalization group invariant 
way 14 , one promotes Gfl avour to a symmetry of the theory by introducing auxiliary fields Yjj and 
Yd transforming under SU(3)g as 

Y v (3,3,1) and Y D (3,1,3). (2.2) 

The Yukawa couplings are then introduced as background fields of these so-called spurions trans- 
forming under the flavour group. An effective theory satisfies the criterion of MFV if all higher- 
dimensional, constructed from SM and Y fields, are invariant under CP and under the flavour group 
Gflavour [l^l - 

In the construction of the effective field theory, operators with arbitrary powers of the dimen- 
sionless Yum have to be considered in principle. However, the specific structure of the SM, with its 
hierarchy of CKM matrix elements and quark masses, drastically reduces the number of numerically 
relevant operators. For example, it can be shown that in MFV models with one Higgs doublet, all 
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FCNC processes with external cf-type quarks are governed by the following combination of spurions 
due to the dominance of the top Yukawa coupling y t : 



(YuYfoij « rfV£V 3j , (2.3) 

where a basis is used in which the d-type quark Yukawa is diagonal. 

There are two strict predictions in this general class of models which have to be tested. First, the 
MFV hypothesis implies the usual CKM relations between b — > s, b —> d, and s —> d transitions. For 
example, this relation allows for upper bounds on NP effects in BR(£? — > X^j), and BR(i? —> X s vv) 
using experimental data or bounds from BR(1? — > X s ^), and BR(if — > tt + vv), respectively. This 
emphasizes the need for high-precision measurements of b — > s/d , but also of s — > d transitions such 
as the rare kaon decay K — > 'nvv. 

The second prediction is that the CKM phase is the only source of CP violation. This implies 
that any phase measurement as in B — > <f)K s or AMb (s/j) is not sensitive to new physics. This is an 
additional assumption because the breaking of the flavour group and the discrete CP symmetry is in 
principle not connected at all. For example there is also a renormalization-group invariant extension 
of the MFV concept allowing for flavour- blind phases as was shown in Ref . [45] ; however these lead to 
non-trivial CP effects, which get strongly constrained by flavour-diagonal observables such as electric 
dipole moments [45] . So within the model-independent effective field theory approach of MFV we keep 
the minimality condition regarding CP. But in specific models like MSSM the discussion of additional 
CP phases within the MFV framework makes sense and can also allow for a natural solution of the 
well-known supersymmetric CP problem, see for example Refs. |23[|24| . 

Scenarios with two Higgs doublets with large tan/3 = 0(mt/mb) allow for the unification of top 
and bottom Yukawa couplings as predicted in grand-unified models and for sizable new effects in 
helicity-suppressed decay models. There are more general MFV relations existing in this scenario due 
to the dominant role of scalar operators. However, since tan (3 is large, there is a new combination of 
spurions numerically relevant in the construction of higher-order MFV effective operators, namely 

(YdY^j « yX , (2.4) 

which invalidates the general MFV relation between b —> s/d and s — > d transitions. 

For more details we refer to the very recent complete mini-review on MFV [32]. Here we only 
add two issues on the application of the MFV hypothesis to the minimal supersymmetric standard 
model (MSSM). Most interestingly, the MFV hypothesis can serve as a substitute for R-parity in the 
MSSM |33,34|. MFV is sufficient to forbid a too fast proton decay because when the MFV hypothesis 
is applied to R-parity violating terms, the spurion expansion leads to a suppression by neutrino masses 
and light-charged fermion masses, in this sense MFV within the MSSM can be regarded a natural 
theory for R-parity violation. Secondly, the MFV framework is renormalization-group invariant by 
construction, however, it is not clear that the hierarchy between the spurion terms is preserved when 
running down from the high scale to the low electroweak scale. Without this conservation of hierarchy, 



the MFV hypothesis would lose its practicability. However, as explicitly shown in Refs. 30,31 , a 
MFV-compatible change of the boundary conditions at the high scale has barely any influence on 
the low-scale spectrum. Finally, the MFV hypothesis solves the NP flavour problem only formally. 
One still has to find explicit dynamical structures to realize the MFV hypothesis like gauge-mediated 
supersymmetric theories. And of course the MFV hypothesis is not a theory of flavour; it does not 
explain the hierarchical structure of the CKM matrix and the large mass splittings of the SM fermions. 



4 



2.2 Effective hamiltonian 

Our analysis is based on the following MFV effective hamiltonian relevant to b — > s transitions (and 
also for b — > d transitions with obvious replacements) |14| : 



^Y.^ y us y ub+v: s v cb )ct + v: s v tb ct]p l+ v: s v tb cx + h. c . (2.5) 



with 
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P9 = 
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P10 - 


^(^ME^'W) • 



(2.6) 



In addition we have the following scalar-density operator with right-handed 6-quark{^] 



e 2 



P^ = Y^(sLbR)(£ R e L ) , (2.7) 



There is no reason to update the MFV analysis of precision s — > d transitions of Ref. [15]. For 
completeness we state that for the rare decays K — > tti>v, we have the following simple effective 
hamiltonian, 

l=e,[i,T 

with the operator P V y and the corresponding possible NP contribution 5C u p, 

Pvv = (S7a.<*) (^7 m (1 - 7s) v t ) , y u = -^- (\ t (X t + bC vV ) + ReX c P UtC ) , (2.9) 

I 'US I ^ ' 



36 



38 



and X q = V* s V qd , X t = 1.464 ± 0.041, P U)C = (0.2248) 4 P U)C U35J and P n , c = 0.41 ± 0.04 

We follow here the analysis of Ref. jl5| and consider NP in the FCNC operators P7,P8,Pg,Pio 
and in the two scalar operators P V y and Pq only; as argued, in principle most of the possible NP 
contributions to the four-quark operators Pi_6 could be reabsorbed into the Wilson coefficient of the 
FCNC operators. The NP contributions to the Wilson coefficients are parameterized as: 



5C t ^ b ) = Cr y (^)-Ct M ^ b ). (2.10) 
where the Cf M ([i b ) are given in Table [l] 



Within the MFV framework, the corresponding Wilson coefficient C is related to Cq 1 = m^Cs and Cq 2 = mjCp 
by Co = 2Cq 1 — —2Cq 2 , where the operators are defined as Qi = Qs/m b = (sl&r)(11), Q2 = Qp/m b — 
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C5(j*) 


-0.2974 


-0.1614 


4.2297 


-4.2068 






Table 1: SM Wilson coefficients at = m^ oe and /xq = 2Mty to NNLO accuracy in a. 



3 Observables and theoretical uncertainties 

We present the various AF = 1 observables which we use in our MFV fit or which we want to 
constrain or predict. We focus on their dependence on the (non-standard) Wilson coefficients of the 
MFV effective theory and on the main sources of the theoretical uncertainty. 



3.1 Radiative decay B — > X Stl ff 



The branching fraction for B — > X q j (q = s, d) for a photon energy cut E^ > Eq can be parameterized 



as 



BR(B -> X ql ) Ei>Eo = BR(S -> X c eu] 



6a e 



exp 



ttC 



Vt q y tb 



V, 



d, 



P(E ) + N(E ) 



(3.11) 



where a P 



a 



on shell 



39 , C = \V, 



ub 



cb\ 



T[B -> X c ev]/T[B -> X u ev\ and P( J E ) and iV(Eb) 
denote the perturbative and nonperturbative contributions, respectively. The latter are normalized 
to the charmless semileptonic rate to separate the charm dependence. The perturbative part of the 
branching ratio of B — > X s ^ is known to NNLL precision 40 , while the nonperturbative corrections 
are now estimated to be well below 10% [42]. The overall uncertainty consists of nonperturbative 



(5%), parametric (3%), perturbative (scale) (3%) and m c -interpolation ambiguity (3%), which are 
added in quadrature. An additional scheme dependence in the determination of the pre-factor C 
has been found [H] ; it is within the perturbative uncertainty of 3% [43]. The dependence of the 
dominating perturbative part from the Wilson coefficients can be parametrized [44] at NNL: 



P(£o) = P (0) (w,) + 



P! (1) (/x 6 ) + P^(E ,fi b ) + O {a 2 M) , where (3.12) 



>(i), 



^ (0) (^) 



a 



(0)eff 



2C? )eff ( w )r- 



(l)eff 
7 



(3.13) 



The functions k!>}' can be found in Ref. 



44 



The effective Wilson coefficients are given in the 
Appendix. We stress that we have used NNLL precision (means inclusion of O (a^ (//&)) terms) in our 
numerical analysis. 

The branching ratio of B — > X^ is only known to NLL QCD precision |45| . The error at this order 
is dominated by a large scale renormalization uncertainty of more than 12% and by uncertainties due 
to CKM matrix elements of 10%. However, in view of the large experimental error the NLL precision 
is still appropriate within our analysis. 
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3.2 Isospin asymmetry Aq(B — > K*j) 



Another important observable which is already measured is the isospin breaking ratio. It arises when 
the photon is emitted from the spectator quark: 



A ± 



(3.14) 



where the partial decay rates are CP-averaged. In the SM spectator-dependent effects enter only 
at the order A /nib, whereas isospin-breaking in the form factors is expected to be a negligible ef- 
fect. Therefore, the SM prediction is as small as 0(8%). Moreover, a part of the A/m^ (leading) 
contribution cannot be calculated within the QCDf approach what leads to a large uncertainty [46| . 
However, the ratio is shown to be especially sensitive to NP effects in the penguin sector. The isospin 



asymmetry can be written as 46 



A = Re(b d - b u ) , 
where the spectator dependent coefficients b q take the form: |^] 



12tt 2 / b Q q 



f ^K 1 + 



fK*rn K * 



K 



2q 



Here the coefficient a'j reads |48 



4tt 



C 2 (fib)G 2 (x cb ) + Cs(fJ-b)Gi 



(3.15) 



(3.16) 



(3.17) 



4tt 



C 2 (fi h )H 2 (x cb ) + C 8 (fi h )H s 



where fi^ = V^W^fe i s the spectator scale. The functions G 2 , G 8 , H 2 , and H$ can be found in Ref. 48) . 
The functions K\ and K 2q can be written in function of the Wilson coefficients Cj at scale fib |46|: 



C 6 (fi b ) + 



C 5 (Mfc 
N 



F± + 



Cf a s (Hb) \ ( m h \' 
N Air I \m B J 



C 8 (fi b )X A 



(3.18) 



-C 2 (fib 



Un r ^ + hF ± -G ± (x cb ) 

6 fib o 



+ T\ 



2g 



VusVub 



v* s v cb 

Cf a s (fib 



C 2 {fi b ) + 



Ci(fib) 



N 



Squ + C 4 (fl b ) + 



CM 

N 



(3.19) 



+ 



N 4tt 



C 2 (fib 



4 . m b 
- In — 
3 fib 



H±(xcb) + r 2 



where x c b 



2 

ml 



3 and Cf = 4/3 are colour factors. The convolution integrals of the 



-f and N 

hard-scattering kernels with the meson distribution amplitudes F± , Gj_ , H± , and Aj_ can be found in 



Ref. 46 , also the residual NLO corrections r\ and r 2 . 



In this subsection we use the Wilson coefficients of the traditional basis 



47 
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Lattice QCD Group Ref. 




Sb 


ETMC-11 


65 






232 


± 


10 MeV 


195 


± 


12 MeV 


Fermilab-MILC-11 


66 


67 




242 




9.5 MeV 


197 


± 


9 MeV 


HPQCD-12 


68 






227 


± 


10 MeV 


191 


± 


9 MeV 


Average 


234 


± 


10 MeV 


194 


± 


10 MeV 



Table 2: Average of lattice QCD results used in this work. 



3.3 Leptonic decays B S ( i — > 

The rare decay B s — > fj,^ '/z~ proceeds via Z° penguin and box diagrams in the SM. It is highly 
helicity-suppressed. However, for large values of tan f3 this decay can receive large contributions. In 
general, within MFV the pure leptonic decay B s — > £ + £~ receive contributions only from the effective 
operators Pio and Pq. These are free from the contamination of four-quark operators, which makes 
the generalization to the b — >■ d case straightforward. The branching fraction in the MFV framework 
is given by 



BR(B 8 



of 



64vr 3 



VtbK 




4m^ 



in 



(3.20) 



4m£ 



mB a 



m b + m s 



K/2) 



+ 



™B a 



m b + m s 



m,. 



™B a 



where f Ba is the B s decay constant, m Ba is the B s meson mass and tb 3 is the B s mean life. 

The main theoretical uncertainty comes from the B s decay constant fs s , which has recently been 
re-evaluated by independent lattice QCD groups of Table [2} Their 4.3% uncertainties agree, as do 
their results within these uncertainties, so that following Ref. |86| we have chosen an average of these 
three results in what follows. This implies a 8.7% uncertainty on the branching ratio. The most 
important parametric uncertainty comes from the CKM matrix element Vt s with 5%. 

Within the MFV scenario the Bj — > £ + £~ rate can be obtained from the one of B s — > £ + £~ with 
the exchange (Vt s , mB 3 , Tts> fs a ) — > (Vtd, m B d , fs d )- This implies a very important MFV relation 
{0(riid/m s ) are neglected), 



fB a m Ba 

T(B d ^£+£-) ~ f Bd m Bd 



is 



td 



(3.21) 



3.4 Inclusive B — > X s ji + n and B — > X s t + t 

The decay B — > X a £ + £~ is particularly attractive because it offers several kinematic observables. The 
angular decomposition of the decay rate provides three independent observables, Ht, Ha and Hl, 
from which one can extract the short-distance electroweak Wilson coefficients that test for NP: 



d 3 T 
dq 2 dz 



3 r 



(1 + z 2 )H T {q 2 ) + 2(1 - z 2 )H L (q 2 ) + 2zH A {q 2 



(3.22) 



Here z = cos0£, 0£ is the angle between the negatively charged lepton and the B meson in the center- 
of-mass frame of the dilepton system, and q 2 is the dilepton mass squared. Ha is equivalent to the 
forward-backward asymmetry, and the dilepton-mass spectrum is given by Ht + Hl. The observables 
mainly constrain the Wilson coefficients Cf , Cg and C£g . 
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One defines perturbatively dominated (means theoretically clean) observables within two dilepton- 
mass windows avoiding the region with the cc resonances: the low-q 2 region (1 GeV 2 < q 2 < 6 GeV 2 ) 
and the high-q 2 region (q 2 > 14 .4 GeV 2 ). 

In order to show the dependence of the observables on the Wilson coefficients, we use the con- 



ventions of Ref. 49 . We note that for the numerical evaluation we have used the conventions in 
Ref. 50 , in particular we have chosen the charmless semileptonic decay rate as normalization. For 
the branching ratio within the MFV framework we find: 

dBR(B -> X s l+t 



\mn >x c iv) 

as A7r 2 f(z)K(z) \V cb 




x \C§°» | 2 (1 + 2 ^)(l + 28) (l + ^r 99 (s)) + 4|Cn 2 (l + 2 ^)(1 + ?) (l + °^t 77 (s) 



2 111 / Cv \ 2 Til / CS> 

\C^ W \ 2 [{1 + 2s) + -^(1 - 4s)] (1 + -r 99 (s)) + \2Re{Cr w Cr w *){l + —^)(l + -^t 79 (s) 



7T / S V 7T 



+ ^| 2 (J - 2m 2 ) - 3Re(C?rC e *)rne + 6%^ + + + 8&% + ^ /d| , (3.23) 



where the hat indicates a normalization by mj. The functions Tj correspond to specific bremsstrahlung 
terms. As indicated, further (but finite) bremsstrahlung, electromagnetic and power corrections have 



to be added, see Ref. [49] for more details. Our formula (3.23) is consistent with the results in Ref. |51|. 



For the dependence of the forward-backward asymmetry on the Wilson coefficients we find within 
the MFV setting: 



A 



fb(s 



* (PER 

l dsdz 



_i dsdz 



-B(B -)• X c lv) 



3a 2 



\v tb v t 



* |2 
ts\ 



4n 2 f(z)K(z) \V cb \ 



1 



! (1 



4m? 



Rc(C^ ew C^ w *)s 



l + -T 9W (s))+2Re{Cr w C n 10 

7T 



i + — mo(s) 

7T 



(3.24) 



+ Re( (C^ ew /2 + C? ew 



Ci*)rh t 



brems I g\ i rem 



(s). 



We note that Eq. (3.24) is consistent with the results in Ref. [51], but disagrees with Refs. [52, 53 



for the scalar contributions. The new Wilson coefficients are defined in Ref. [49] and are given in the 
Appendix. 

In the low-q 12 region the theoretical uncertainty is around 7% for the branching ratio, however 
there is an additional 5% uncertainty due to nonlocal power corrections to be added [50]. In the high- 
q 2 region, one encounters the breakdown of the heavy-mass expansion at the endpoint. However, for 
an integrated high-q 2 spectrum an effective expansion exists in inverse powers of m^ ff = m& x (1 — 
y/s m i n ) rather than ra b . The resulting large theoretical uncertainties in the high-q 2 due to the power 
corrections of around 25% could be significantly reduced by normalizing the B — > X s £ + £~ decay rate 
to the semileptonic B — > X u iv decay rate with the same q 2 cut [54]. For example, the uncertainty 
due to the dominating l/m| term would be reduced from 19% to 9% [50| . 
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3.5 Exclusive decay B — > K*( 



The exclusive semi-leptonic penguin modes offer a larger variety of experimentally accessible observ- 
ables than do the inclusive ones, but the hadronic uncertainties in the theoretical predictions are in 
general larger. 

The physics opportunities of B — > K*U (£ = e, fx, r) decays depend strongly on the measure- 
ment of their angular distributions. This decay with K* on the mass shell has a 4-fold differential 
distribution [57[[58] 



d A Y[B -»• K*(-> Ktt)U] _ 9 
dq 2 d cos 9i d cos 6 k dip 32ir 



^2j t (q 2 )gi(9 h e K ,, 



(3.25) 



w.r.t. the dilepton invariant mass q and the angles 6i, 6k, and <p (as defined in [55]). It offers 12 
observables Ji(q 2 ), from which all other known ones can be derived upon integration over appropriate 
combinations of angles. 

The J{ depend on products of the eight theoretical complex K* spin amplitudes A{, A^^ , A t , As- 
The Ji are bi-linear functions of the spin amplitudes such as 



J 



\A^\ 2 + \Ak\ 2 + \A?\ 2 + \A$\ 2 



(3.26) 



with the expression for the eleven other Jj terms given for example in 56 , 58 59 

The dilepton invariant mass spectrum for B — > K*£ + £~ can be recovered after integrating the 
4-differential distribution over all angles, while the (normalized) forward-backward asymmetry ^4fb 
can be defined after full <p and 9k* integration [60] (Ji = 2Jf + Jf): 



dq 2 



■h 

3 



^4fb(c? 2 ) 



dcost 



d 2 T 



dT 



dq 2 d cos 6i / dq 2 



3 J 6 
8dT/dq 2 



(3.27) 



Moreover, the fraction of the longitudinal polarized K* is given by Fl = (3 Jf — J|)/(4dr/cig |2 ).These 
three observables represent the early ones, which have been measured already by the B factories and 
now with much better precision by the LHCb experiment. 

With more luminosity, theoretically much cleaner angular observables will be available. In the 
low- and high-g 2 region it is always appropriate to design optimized observables by using specifically 
chosen normalizations for the independent set of observables. In the low-g 2 region, specific ratios of 
observables allow for a complete cancellation of the hadronic uncertainties due to the form factors in 



leading order and, thus, for a high increase in the sensitivity to new physics structures [55,59 , for 
example the transversity amplitudes: 



.4 



(2) _ 1^3 

" 2 J| 



(3)_ / 4J 2 + /? 2 J 2 
V -2J 2 c (2Jf + J 3 



A 



A 



(4) 

T 



iJi + PfJ 2 ' 



(3.28) 



In the high-g 2 region, two groups of ratios of o bservab les can be constructed which dominantly depend 
either on short- or on long-distance physics 



transversity observables were proposed: 
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(i) 

In addition to the A T observables some new 



H. 



(i) 



V2J 4 



y/-j°(2jz-j 3 y 



H. 



(2) 



Plj5 



y/-2 J| (2 J| + J 3 ) 



H, 



(3) 

T 



PiJs 



4 



(3.29) 
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In the high-g 2 region H^'^ depend only on short-distance information in leading order, while Fl and 



.4 



(2,3) 



, depend only on long-distance quantities. ^ 
The theoretical treatment in the low- and high-g 2 is based on different theoretical concepts. Thus, 
the consistency of the consequences out of the two sets of measurements will allow for an important 
crosscheck. 

In the \ow-q 2 region, the up-to-date description of exclusive heavy-to-light B — > K*£ + £~ decays 
is the method of QCD-improved Factorization (QCDF) and its field-theoretical formulation of Soft- 
Collinear Effective Theory (SCET). In the combined limit of a heavy 6-quark and of an energetic 
K* meson, the decay amplitude factorizes to leading order in A/ nib an d to all orders in a s into 
process-independent non-perturbative quantities like B — > K* form factors and light-cone distribution 
amplitudes (LCDAs) of the heavy (light) mesons and perturbatively calculable quantities, which are 
known to O(al) 16 17 . Further, the seven a priori independent B — > K* QCD form factors reduce 
to two universal soft form factors £_ui [64]. The factorization formula applies well in the range of the 
dilepton mass range, 1 GeV 2 < q 2 < 6 GeV 2 . 

Taking into account all these simplifications the various K* spin amplitudes at leading order in 
-^-QCD/^ft and a s turn out to be linear in the soft form factors and also in the short-distance 
Wilson coefficients which allows to design a set of optimized observables in which any soft form factor 
dependence (and its corresponding uncertainty) cancels out for all low dilepton masses q 2 at leading 
order in a s and Aqcd/ 77 ^ 



55 
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A^ R = V2Nm B (l - s) 



(4 eS) =F do) + 



2m& (eff) 

S L 7 



K* 



(3.30a) 



-V2Nm B (l - s) 



(cf ] Tc 10 ) + 2 -^cr 



(3.30b) 
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L,R 



NrriB 



:(1 



S? 



(4 eff) =f Cio) + 2m b C 



(off) 



A'* 



(3.30c) 



Nnis 



(1 



-10 



-a 



£|l (Ek* 



(3.30d) 



A, 



Nm 2 B ( 
2niK*Tnb 



■i)cg 



€\\(Ek*) , 



(3.30e) 



with s = q 2 /m 



Bi 



111; 



mi/niB- Here we neglect terms of 0(fn 2 K ,\ 



numerical analysis. The factor collects all pre-factors and can be found in Ref. 1 55 59 



but we include these terms in our 

The soft 



form factors are fixed in a specific factorization scheme using QCD sum rule techniques as discussed 
in the Appendix. 

However, in the early observables, namely dT / dq 2 , Ap~B , Fl there is still a large theoretical uncer- 
tainty due to the form factors which do not cancel out to first order in these cases. 

Within the QCDF/SCET approach, a general, quantitative method to estimate the important 
Aqcd /mb corrections to the heavy quark limit is missing. In semileptonic decays a simple dimensional 
estimate of 10% is often used. Under the assumption that the main part of the AQco/^b corrections 
is included in the full form factors, the difference of the theoretical results using the full QCD form 



There are three observables which are already measured beyond the early observables mentioned above: S3 = 
(J3 + J3)/[d(r + f )/dq 2 ], Aim, and the isospin asymmetry, but all three observables have no significant impact on the 
MFV scenario yet. 
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factors on one hand and the soft form factors on the other hand confirms this simple dimensional 
estimate. In fact, the comparison of the approaches leads to a 7% shift of the central value. 

The low-hadronic recoil region is characterized by large values of the dilepton invariant mass 

It is shown that local 
K*£+£~ 



q 2 > (14 — 15) GeV 2 above the two narrow resonances of J/ip and ip(2S) 
operator product expansion is applicable (q 2 



~ m b) 



19 



and it allows to obtain the B 



matrix element in a systematic expansion in a s and in A/m&. Most important, the leading power 
corrections are shown to be suppressed by (Aqcd/^&) 2 or « s Aqcd /mb 19 and to contribute only 



at the few percent level. The only caveat is that heavy-to-light form factors are known only from 
extrapolations from LCSR calculations at low-g 2 at present. But this may improve in the future when 
direct lattice calculations in the high-t/ 2 are available 
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There are improved Isgur-Wise relations between the form factors in leading power of A/mj,. Their 
application and the introduction of specific modified Wilson coefficients lead to simple expressions for 
the K* spin amplitudes to leading order in l/m^ in the low recoil region, for example we have pi] 



A 



L,R 



+i 



/-^effjinod fi 



) + K— C 7 fl 



(3.31) 



> <{C 9 =f c w) + K — c 7 I /||> 



A 



L,R 



a 



eff ,mod 
9 



TC 10 )+k 



-»*6 £,cff,mod j 



/0) 



(3.32) 



(3.33) 



where the form factors f±, fn, and /o are linearly connected to the QCD form factors (see Ref. |61|). 
The modified effective Wilson coefficients introduced in Ref. |18| are given in the Appendix. Then, 
the three considered observables at leading order can be written in the high-g 2 region as 
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dr 

dq 2 



2pi x (/ 2 + /|+/f), Afb = 3 



Pl x M F] 

Pi ui + fl +/,?)' 



/o 2 



/ 2 +/ 2 +/,r 



(3.34) 



'|| ^ Jo 

where only the two independent combinations of Wilson coefficients enter, namely 

2 



Pi 



2rhb 



Cg -|- K - C-j 



cff 



+ |C 10 r , P2 = Re{ [ Cf + K^Cf ) C 10 



(3.35) 



pi and p2 are shown to be largely /U-scale independent 18 



As mentioned above, the leading power corrections of the OPE arise at 0(a s A/mi,,m^/Q ) and of 
the order of a few percent. The A/mj, corrections to the amplitudes from the form factor relations are 
parametrically suppressed as well, by small dipole coefficients, such that one can estimate the leading 
power correction from the form factor relations to the decay amplitudes as order (2C^ S /Cg S )A/mb- 
So in general, the dominant power corrections to the transversity amplitudes are of the order of a few 
percent |61|. 



4 Numerical details 

Within our numerical analysis we use the most recent LHCb results for the exclusive decays B s — > 
and B — > K*p + p~, and Belle, Babar and CDF results for the other decays. The experimental 
values are provided in Table [3j For comparison, we also consider the pre-LHCb data as given in 
Table H 
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Observable 


Experiment 


SM prediction 


BR(B -> X s7 ) 


(3.55 ± 0.24 ± 0.09) x 10 -4 77 




(3.08 ±0.24) x 10" 4 


A (B ->■ K*7) 


(5.2 ± 2.6 ± 0.09) x 10" 2 77 


(8.0 ±3.9) x 10 -2 


BR(B -j. X d7 ) 


(1.41 ± 0.57) xj_0" 5 


78 


79 


(1.49 ±0.30) x 10 _i> 


BR(_B S -> n + 


< 4.5 x 10 -9 80] 


(3.53 ±0.38) x 10 _a 


(dBR/dq\B -> ^> + /0} 9 2 eIli6]GeV 2 


(0.42 ±0.04 ±0.04) x 10"' 


81 




(0.47 ±0.27) x 10" ' 


(dBR/dq 2 (B -¥ K* ji + ^")}<j2 6[1 4. 18il6]GeV 2 


(0.59 ±0.07 ±0.04) x 10"' 


81 




(0.71 ±0.18) x 10"'' 


(Afb(B -> #*A« + /0)« 2 e[i,6]GeV 2 


-0.18 ±0.06 ±0.02 181 




-0.06 ±0.05 


(Afb(B -> -^V + M _ )) 9 2 e[14.18,16]GeV2 


0.49 ±0.06 ±0.05 81 


0.44 ±0.10 


gg(AFs(S-»>-K:*At + M _ )) 


4.9±i; L 3 GeV 2 81 


4.26 ± 0.34 GeV 2 


-» #"W~)),J 2 €[l,61GeV 2 


0.66 ±0.06 ±0.04 [8~l]_ 


0.72 ±0.13 


BR(_B -> X s M + M~) 9 2 e [i i6 ] G eV 2 


(1.60 ±0.68) x 10 _7r " 


82 


83 




(1.78 ±0.16) x lO"" 


BR(B — ► X s (i + n ) 9 2 >14 .4Gcv 2 


(4.18 ± 1.35) x 10"' 


82 


83 




(2.19 ±0.44) x 10" ' 



Table 3: Input observables: The experimental data represent the most recent one. The updated SM 
predictions are based on the input parameters given in Table [5] 



There is a remark in order regarding the branching fraction of B s —¥ Its value provided by 

the experiments corresponds to an untagged value, while the theoretical predictions are CP-averaged. 

, the untagged branching ratio is related to the CP-averaged one 



As pointed out recently in 87 
by: 



where 



1 + ^Ar Vs 



1 



vi 



: t b AT s = 0.088 ±0.014 



(4.36) 



(4.37) 



and 



with 



Aai 



S = 



|P| 2 cos(2</> P ) - \S\ 2 cos(2<p s ) 



\P\ 2 + \S\< 



mf, 

4 



Ml 



C l /2 



Ml 



2m M m b 



Cio 
°io 



M 



+ 



1 



m s Cf M 



-C&/2 



2mn mi, + m s 



°io 



and 



ip s = arg(S) , (f P = arg(P) . 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



The obtained branching ratio can then be directly compared to the experimental result. 

The SM predictions entering the MFV-fit are based on the theoretical analyses given in Section [3] 
We have used the input parameters of Tableland the program Superlso v3.3 [75,76 in order to 
update SM predictions. They are given in Table [3] 

To obtain constraints on the Wilson coefficients, we scan over 6C7, SCr, 8Cq, <5Cm and SCq. For 



each point, we then compute the flavour observables using Superlso v3.3 75,76 and compare with 
the experimental results by calculating x 2 as: 



X 



E 



(or - q\ 

K exp ) 2 + (* 



th\2 



th\2 



(4.42) 
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Observable 


Experiment 


BR(B -> X 3 >y) 


(3.55 ± 0.24 ± 0.09) x 10" 4 77 




(5.2 ±2.6 ±0.09) x 10" 2 77 


BR(B ->■ X d j) 


(1.41 ±0.57) x 10" b 78 


79 


BR(B S ->■ fi+fi-) 


< 5.8 x 10~ 8 84 


{dBR/dq'(B ^*£ + ^-)} 92e[1 , 6]GoV 2 


(0.32 ±0.11 ±0.03) x 10 _Y 


85 




{dBR/dq J (B -J- ■K'*M + M~)} l; 2e[i4.i8,i6]Gev2 


(0.83 ±0.20 ±0.07) x 10 _Y 


85 




(,4fs(-B -> -R"*M + M _ )) 9 2 e [i i6 ] G ov2 


0.43 ±0.36 ±0.06 85 


(,4fs(-B -> ^*M + M~))q 2 S[14.18,16]GcV 2 


0.42 ±0.16 ±0.09 || 


{Fl(B -¥ ^*M + M~)) 9 2 e[i,6]Gev 2 


0.50 ±0.30 ±0.03 85 


BR(B -> X s /i + M"),2 e[16]GcV 2 


(1.60 ±0.68) x 10~ D |82 


83 




BR(B — > X s fi + fi~) q 2 >14AGcV 2 


(4.18 ± 1.35) x 10" Y p2 


83 





Table 4: Experimental data used for pre-LHCb fit. 



m B = 5.27950 GeV 
m K * = 0.89594 GeV 


71 
71 




m Bs = 5.3663 GeV 
|^*|=0.0403i°;°^ 


71 
71 




mf s (m 6 )= 4.19l° 1 £ 8 i GeV 

m Vole = g ± q g ± q g 


71 
71 




mf s (m c ) = 1.291° ?! GeV 
m M = 0.105658 GeV 


71 
71 




cts(Mz) = 0.1184 ±0.0007 

a s (p b ) = 0.2161 

sin 2 6 w (Mz) = 0.23116(13) 


71 
71 




&em{M Z ) = 1/127.916 
&em(m b ) = 1/133 

G F /{hcf = 1.16637(1) GeV" 2 


71 
71 




f B = 194 ± 10 MeV Tabic I 
f Bs = 234 ± lOMeV Table \ 


> 
I 


t b = 1.519 ± 0.007 ps m 
t Bs = 1-472 ± 0.026 ps WT 




/k*,±(1 GeV)= 0.185 ± 0.009 GeV 
a li± (l GeV)= 0.10 ±0.07 
a 2 ,x(l GeV)= 0.13 ±0.08 
A B ,+(1 GeV)= 0.46 ± 0.11 GeV 


72 
73 
73 
74 




= 0.220 ±0.005 GeV 
01,11(1 GeV)= 0.10 ±0.07 
a 2 ,ll(l GeV)= 0.09 ±0.05 


72 
73 
73 




pole 

fi b = ml 

Hj = VO-5 x Mb GeV [l7] 


Mo = 2Mw 



Table 5: Input parameters. 



where 0^ xp and O* are the central values of the experimental result and theoretical prediction of 
observable i respectively, and cr^ xp and af 1 are the experimental and theoretical errors respectively. 
The global fits are obtained by minimization of the \ 2 - For LHCb fits, we consider the measurements 
of the observables given in Table[3]while for the pre-LHCb the measurements of Table|4]are considered. 



5 Results 

5.1 Separate bounds 

We first study the individual constraints from the observables described in Section [3| The main 
players in our analysis are the radiative decay B — > X s ^, the leptonic decay B — > fj- + fi~, and the 
semileptonic decays B — > X s /K* fj,~ /i~ 

Figure [l] shows that similar zones are probed by the inclusive decays B — > X s j and B — > X^. 
The bounds in the (<5CV, 5C%) planes induced by the two inclusive decays are nicely consistent with 
each other as expected in the MFV framework which predicts a strong correlation between the two 
decays. Clearly, due to the smaller theoretical and experimental error the B — > X s j bound is much 
stronger. 
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In the previous MFV analysis [15] the approximation was used that the NP contributions to 
the electromagnetic and chromomagnetic operators appear in a fixed linear combination, namely 
5C-j + 0.35Cs- This additional assumption was necessary in the previous analysis due to the limited 
number of independent experimental measurements. The correlations between 5Cj and 6C$, shown 
in Figure [TJ do not support this simplifying assumption. 




I I 1,1, c I I , 

.5 0.5 1 -U5 0.5 

5 C 7 8 C 7 



Figure 1: 68% and 95% CL bounds on SC-j and 5C% induced by the inclusive decays B — > X s j (left) 
and B — > (right). 



The isospin asymmetry in the exclusive mode B —> K*j brings complementary information to the 
inclusive branching ratios. Figure [2] shows that the isospin asymmetry seems to favor opposite signs 
for bC-j and 5C$. 




Figure 2: 68% and 95% bounds on on bC-j and bC% induced by the isospin asymmetry in B — > K*j. 



The leptonic decays B s — > fi + and Bj — > fJ- + n~ are sensitive for <5Cio and the scalar contribution 
5Cq. The shapes in the corresponding correlation plots induced by the two leptonic decays are very 
similar, thus, highly consistent with each other as can be seen in Figure [3j This feature strongly 
supports the MFV hypothesis which predicts a strong correlation between these two decays as given 
in Eq. (3.21). Of course, the experimental limit for the decay B s — > fi^ is much tighter and 
therefore the present constraints are much stronger. We therefore take the decay B^ —> out of 

the global MFV fit and will make a prediction for this decay within the MFV framework below. 

We notice that the constraint on the scalar coefficient induced by the decay B s — > fi + fi~ is very 
strong and a large scalar contribution is not allowed anymore. 
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Figure 3: 68% and 95% CL bounds on on 5C\g and 5Cq induced by the decays B s — > fi + fi (left), 
B d -)■ n + fi- (right). 



The low-q 2 data of the inclusive decay B — > X s ^ + fi~ and of the exclusive decay B — > K*[i + [i~ 
have similar constraining power, as can be seen in Figure[4j It is nontrivial that the correlation plots of 
the various 5Ci look almost identical for the inclusive and the exclusive mode. B — > X s /j, + has small 
theoretical errors and large experimental errors, while the situation is reversed for B — > K* [i + /i~ . A 
statistical combination of both allows to enhance their effect. However, one realizes the potential of 
the inclusive mode if one takes into account the fact that the recent Babar and Belle measurements 



of the inclusive branching ratios 82 , 83 only use less than a quarter of the available data sets of the B 
factories. The constraints on C\q are similar to those from B s — > fj> + ^ , but contrary to B s —> fi + fi~, 
the constraints on the scalar contributions here is very weak. 

Finally, we note that the allowed values of 5Cg and 5C\o are much smaller in specific NP models 
than within a model-independent analysis, so for example the structure of the CMSSM already bounds 
their values significantly before any experimental data is used (see Ref. |86|). 



5.2 Fit results 

We made two global MFV fits in order to make the significance of the latest LHCb data manifest, 
see Figure [5} First we have used the experimental data before the start of the LHCb experiment 
(pre-LHCb, right plots). These measurements are listed in Table |4j Then we have included the latest 
LHCb measurements given in Table [3] (post-LHCb, left plots). 

Here Cs is mostly constrained by B — > X s ^-, while Cj is constrained by many other observables 
as well. Cg is highly constrained by b — > s/j, + fj,~ (inclusive and exclusive). C\q is in addition further 
constrained by B s — > [i + . Cg is dominantly constrained by B s —¥ [i + [iT . 

There are always two allowed regions at 95% CL in the correlation plots within the post-LHCb fit; 
one corresponds to SM-like MFV coefficients and one to coefficients with flipped sign. The allowed 
region with the SM is more favored. The various JCj-correlation plots show the flipped-sign for C-j 
is only possible if Cg and Cio receive large non-standard contributions which finally also change the 
sign of these coefficients. 

We have also studied the impact of the LHCb measurements of the branching ratio, of the forward- 
backward asymmetry, and of the K* polarization within the exclusive decay B — > K*l + i~ by taking 
these LHCb measurements out of the fit. The results in Table [6] show that these pieces of experimental 
information from the LHCb experiment are very important. They significantly reduce the allowed 
areas for 6Cg and 5C\q. 
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Figure 4: 68% and 95% CL bounds on various 5Ci induced by the decays B — > X s fi + fi (left) and 
B -> K*/j, + fi~ (right) at low-g 2 . 



5.3 Predictions within the MFV benchmark 

With the help of the results of the global fit, which restricts the NP contributions <5Cj, we can now 
derive several interesting predictions of observables which are not yet well measured. This analysis 
allows to spot these observables which still allow for relatively large deviations from the SM even in 
the MFV benchmark scenario. 

• For the branching ratio of the decay B — > X s tt we get the 95%CL bounds 

0.2 x 1(T 7 < BR(B -»• V s r + r-) (?2>14 4GcV 2 < 3.7 x 1(T 7 . (5.43) 
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Figure 5: Global MFV fit to the various NP coefficients 8C% in the MFV effective theory with (left) 
and without experimental data of LHCb (right). 
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Figure 6: Global MFV fit with the latest data set excluding all LHCb measurements of B — > K* \x 
observables. 



This has to be compared with the SM prediction 

BR(£ X s t + t~) 2>UAGeY 2 = (1.61 ± 0.40) x 10~ 7 . 



(5.44) 



So there are still large deviations from the SM prediction of this observable possible within the 
MFV scenario. And as stated above, any measurement beyond the MFV bounds would indicate 
the existence of new flavour structures. 

• For the zero-crossing of the forward-backward asymmetry in the inclusive decay B — > X s fj, + fi~ , 
we get the lower bound at the 95%CL 



A FB {ql) =0; 1.94 GeV 2 < q 2 



(5.45) 



while the very precise SM prediction is (<?o) SM = (3.40 ± 0.25)GeV 2 . There is natural upper 
bound given by the cut due to the charm resonances. Due to the theoretical cleanliness of this 
observable, there are still large deviations from the SM prediction possible within the MFV 
benchmark. This is also true for the the complete function Afb(q 2 )- 

We have taken the measurements of the decays B — > and B d — > \i~ out of the global fit 
and find the following MFV predictions again at the 95%CL: 



(5.46) 



1.0 x 10~ 5 < BR(B -> X d -f) < 4.0 x 10~ 5 ; BR(B d < 3.8 x 10~ 10 . 

The corresponding SM predictions are: 
BR(S^V d7 ) SM = (1.49±0.30)xl0~ 5 ; BR(B d -»• n + ^) SM = (1.11±0.27) x 10" 10 , (5.47) 
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and the present experimental data (see Table [3]) is: 

BR(B A d7 ) Exp - = (1.41 ± 0.57) x 10~ 5 ; BR(£ d -> /x + ^) Exp < 10.0 x 10~ 10 . (5.48) 

So the present -B — )• ^7 measurement is already below the MFV bound and is nicely consistent 
with the correlation between the decays B — > X s ~f and B — > predicted in the MFV scenario. 
In the case of the leptonic decay B^ — > fi + fi~, however, the MFV bound is stronger than 
the current experimental limit. And there are still sizable deviations from the SM prediction 
possible within the MFV but an enhancement by orders of magnitudes due to large tan j3 effects 
are already ruled out by the latest measurements. 

For the large set of angular B — > K*n + fj>~ observables discussed in Section [3] we also can easily 
derive their MFV predictions. In an exemplary mode, we give the 95%CL MFV predictions for 
the Ay) , averaged over the low-q 2 region (lGeV 2 < q 2 < 6GeV 2 ), 

- 0.065 < (A^) < -0.022; 0.34 < (A$p) < 0.99; 0.19 < (A^) < 1.27, (5.49) 

and for the , averaged over the high-g 2 region (14.18 GeV 2 < q 2 < 16GeV 2 ), 

(H^) = 1, -1.01 < {H^} < -0.44, -1.01 < (i4 3) > < -0-44. (5.50) 

Due to the experimental and theoretical uncertainties of the A^ observables, discussed in 



Refs. 1 55 1 or [63], the predicted MFV range cannot be really separated from the SM prediction; 
any significant deviation from the SM prediction indicates new flavour structures. But for the 

(i) 

H)p observables deviations form the SM are still possible within the MFV scenario. 

For the rare s — >■ d transitions we refer again to the previous analysis in Ref. |15| . We just 
state that the existing rather weak experimental bound on the branching ratio of BR(A' + — >• 
implies MFV predictions for the flavour observables BR(1? — > K^vv) and BR(A^ — > 



7T ' vv 

r0 



89 



ttvv). Moreover, as in the case of the two decays B s c [ — > [ifi, the charged and neutral K — > ■nvv 
decays are governed by only one parameter, namely the real coefficient 5C u p in the MFV effective 



theory (see Eq.(2.9)), thus, the ratio of the two K — > ttuu allows for an important model- 
independent test of the MFV hypothesis. There will be two dedicated kaon experiments, E-14 
Koto at J-PARC [90] and Na42 at CERN (9l], for this task in the near future. 
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A Effective, new, and modified Wilson coefficients 

The standard effective Wilson coefficients are defined as follows: 

6 6 

Cf(M) = a(M), fori = l,...,6, Cf( f i) = C 7 ( f i) + Y,y j C J (ri, Cf (/i) = tf 8 ( M ) + £ z&Qi). 

3=1 3=1 

(A.l) 
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,20 



then the 



In the MS scheme, one fixes y = (0, 0, — g, — g, — , — ™) and i* = (0,0,1, 
leading-order 6 — )• 57 and b ^ sg matrix elements of the effective Hamiltonian are proportional to the 
leading-order terms in Cy ff and C| ff 44 . Moreover, we have C^q = Cio(/i), and Cg ff is defined as: 



Cf(s) = <7 9 (/z) + ^ dbfo® In ( ^ ) + ^ 3 ( M ) + ^C 5 (/x) + ^0>(m) 



i=l 



?7lft\ 4 



64 



64 



27 



+ g (m c , s) ( -Ci(^) + CM + 6C 3 (/x) + 60C 5 (/x) 



7 



32 



+ a) --C 3 (p) - -C 4 (m) - 38C 5 (/i) - -C 6 (/i) 



12 32 
+ //(U. *) ( —C 3 (m) - -C 4 ( M ) - 8C5M - y C 6 ( M ) 



2 " vr/ 3 
where the function g(rh c , s) is given as 



(A.2) 



/ \ 4 1/N 8 16 z 2 / Az 
g(z,s) = -- H z) + - + ----(2+- 



Az 



x < 



2 arctan 
In 



iz-s 



for s < 4 z . 
i 7r for s > 4z . 



(A.3) 



T/ie new Wilson coefficients are introduced in Ref. |49| and are defined as: 



7T 



-fCff 



47T - 



c[ 0) Fi 7 \s) + C^"(s) + c; ll(U 'f^(s 



,(0) P (7) 



€ff(0) p(7), 



Cf-( S ) = (l + C$. 



47T 



(0) P (9), 



,eff(0) P (9), 



(A.4) 
(A.5) 
(A.6) 



The virtual corrections to 0\p and 0% are embedded in F^ 7 ^ and i ? g 7,9 ' ) . The <7j functions indicate 
certain bremsstrahlung contributions (see Ref. [49| ) . 



Finally, the modified effective Wilson coefficients introduced in Ref. 18 include contributions of 



the four-quark operators but also of the gluon dipole operators. It is important to note that these 
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n 


r 2 m 2 R [GeV 2 ] 


mj it [GeV 2 ] 


V 


0.923 


-0.511 5.32 2 


49.40 


At 




0.290 


40.38 


A 2 


-0.084 


0.342 


52.00 



Table 6: Fit parameters describing the q 2 dependence of the form factors V and A\i in the LCSR 
approach [73] . 



quantities are different from the effective Wilson coefficients Cj\ introduced before. 



a 



eff,mod 



C 9 + h(0, q- 
U(m b ,q>) 



C\ + C2 + ~2 ^3 ~~ 2 C4 + 52 C; 
64 



32 



7C 3 + -C i + 76C 5 + -C 6 



+ 



+ 



4ir 



c + — a + — a 

u 3 1 3 5 9 6 



Ci (S(g 2 ) +4C(g 2 )) -3C 2 {2B(q 2 ) - C{q 2 )) - C 8 F^(q 2 ' 



+ 



^C l + ±C 2 + 2C 3 + 20C 5 



^reff,mod /~< 1 

° 7 " ° 7 " 3 



C 3 + i<7 4 + 20C 5 + yC 6 



The functions A, 5, C and F 8 {7) , F 8 (9) are given in Refs. [92] and 



+ 



47T 



16 



(A.7) 



(C x - 6 C % ) A(q 2 ) - C 8 F 8 (7) (g 2 ) , (A. 



B Determination of the soft form factors 

To obtain the soft B — > K* form factors we have used the following factorization scheme 17 

M? 2 ) = M M * V(q 2 ) , (B.l) 
Mb + m K * 

2 M B + m K * 2 M B - mjf. 2 

) = 2 ^ ) " Mi3 A2(g } • (B - 2) 

The full form factors V and j4i )2 have been taken from light-cone sum rule (LCSR) calculations |73|: 



V(q 2 ) 
Al(q 2 ) 
A 2 (q 2 ) 



+ 



r-2 



1 — q 2 /m 2 R 1 — q 2 /m 2 it ' 

7"2 

1 - q 2 /m 2 ' 



+ 



r-2 



1 - q 2 /m 2 it (1 - q 2 /m 2 it ) 2 ' 



(B.3) 
(B.4) 
(B.5) 



where the fit parameters ri )2 , m 2 R and ?7i 2 j t are given in Table |6j 
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